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We calculated for the nearest active galactic nucleus (AGN), Centaurus A (Cen A), the 
flux of high energy cosmic rays and of accompanying secondary photons and neutrinos 
expected from hadronic interactions in the source. We used as two basic models for the 
generation of ultrahigh energy cosmic rays (UHECR) shock acceleration in the radio 
jet and acceleration in the regular electromagnetic field close to the core of the AGN, 
normalizing the UHECR flux to the observations of the Auger experiment. Here we 
compare the previously obtained photon fluxes with the recent data reported by the 
Fermi LAT and H.E.S.S. collaborations. In the case of the core model, we find good 
agreement both for the predicted spectral shape and the overall normalization between 
our earlier results and the H.E.S.S. observations for a primary proton energy dN/dE oc 
E~°' with a ~ 2 or smaller. A broken-power law with high-energy part a = —2.7 leads to 
photon fluxes in excess of the Fermi measurements. The energy spectrum of the photon 
fluxes obtained by us for the jet scenario is in all cases at variance with the H.E.S.S. and 
Fermi observations. 
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1. Introduction 

Progress in cosmic ray (CR) physici^^ has been hampered for long time by the de- 
flection of charged cosmic rays in magnetic fields, preventing the identification of 
individual sources. This problem could be solved by using the neutral messengers 
that should be produced as secondaries in hadronic CR interactions close to the 
source. However, the secondary photons generated are difficult to disentangle from 
photons produced by synchrotron radiation or inverse Compton scattering of elec- 
trons. Moreover, high energy photons are strongly absorbed both in the source and 
propagating over extragalactic distances. By contrast, the extremely large mean free 
path of neutrinos together with the relatively poor angular resolution of neutrino 
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telescopes and the small expected event numbers makes the identification of ex- 
tragalactic sources challenging using only the neutrino signal. Performing neutrino 
astronomy beyond the establishment of a diffuse neutrino background requires there- 
fore most likely additional input, either timing or angular information from high 
energy photon or CR experiments. 

The recently announced evidencJ^I for a correlation of the arrival directions of 
UHECRs observed by the Pierre Auger Observatory (PAO) with active galactic 
nuclei (AGN) may provide a first test case for successful "multi-messenger astron- 
omy." In particular, Ref. 2 finds two events within the search bin of 3.1° around 
the nearest active galaxy, Centaurus A (Cen A). This FR I radio galaxy is located 
close to the supergalactic plane at a distance of about 4 Mpc (see Ref. [3] for more 
details). At present this correlation, not confirmed by the HiRes experiment^], has 
only 3ct C.L. and other source types that follow the large-scale structure of matter 
would also result in an excess of events along the supergalactic plane. 

In Ref. [7l we studied therefore the possibility to observe Cen A using high-energy 
photons and neutrinos. Taking the correlation signal at face value, we used the PAO 
results as normalization of the CR flux and calculated the fiux of accompanying 
secondary photons and neutrinos expected from hadronic interactions in the source. 
Since both the Fermi LAT collaboration^ and the H.E.S.S. coUaborationl^l reported 
recently the discovery of 7-ray emission from Cen A, we have now the opportunity 
to check our predictions against these observations. 

2. Assumptions 

We calculated the flux of high energy cosmic rays and of accompanying secondary 
photons and neutrinos expected from Cen A for two different scenarios: Accelera- 
tion close to the core, either in accretion shocks or regular electromagnetic fields, 
and acceleration in the radio jet. In the first case UV photons are the most impor- 
tant scattering targets. We modeled the primary photon field around the AGN core 
guided by the simplest possible theoretical model - , namely the thermal emission 
from a geometrically thin, optically thick Keplerian accretion disc. The interaction 
depth for photo-hadron interactions can reach Tp^ ^ few. According to the obser- 
vational data we found that pp interactions of UHE protons with the gas provide 
the main source of CR interactions in the second scenario. In this case, moreover, 
diffusion in the turbulent magnetic fields will increase the interaction depth at low 
and intermediate energies. We considered also three spectra dN/dE oc E~°' of the 
injected protons: Power-laws with a = 1.2 and a = 2, and a broken power-law with 
a = 2.7 for E > Eb = 10^^ eV. 

We based our calculations on several simplifying assumptions like the use of an 
one-dimensional geometry and the omission of the acceleration process. In particu- 
lar, we only postulated that acceleration to 10^" eV is possible in the environment of 
Cen A, without demonstrating it for a concrete model. Finally, hadronic interactions 
are simulated with an extension of the Monte Carlo code described in Ref. 13 
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3. Results versus observational data 

Figure [1] displays the particle fluxes predicteci^ from Cen A as function of the energy, 
assuming that the two events observed by PAO around Cen A indeed originate from 
this AGN. The case of acceleration close to the core is shown on the left, while the 
case of acceleration in the jet is shown on the right. From the top to the bottom, 
spectra are displayed for a broken power-law, a = 2, and a ~ 1.2. In addition 
to the injected proton flux (black solid line), we show the flux of protons (black 
dashed), photons (blue solid) and neutrinos (red solid) arriving on Earth. Note that 
the cutoff in the neutrino and proton spectra below 100 GeV is artificial, since we 
neglect neutrinos and protons with lower energies in our simulation. 

In the core model the final proton flux is reduced by photon-proton interactions 
by a factor w 2 above the threshold energy ~ 10^^ eV (left), while diffusion in 
the jet increases the interaction depth for lower energies (right), resulting in the 
effective production of secondaries. Since the CR spectra are normalized to the 
integral UHECR flux above E'th = 5.6 x 10^^ eV, steeper spectra result in larger 
secondary fluxes at low energies. 

While we could show ear lieJ^ these fluxes only together with an upper limit from 
H.E.S.S. and the estimated sensitivity of Fermi for point sources, we have updated 
now these figures with the recently published results from FermP' and H.E.S.S. 
citeAharonian:2009xn. Remarkably, the photon flux in the Fermi and the H.E.S.S. 
energy range has approximately the same power-law exponent {a ^ 2.6 and 2.7), 
but the latter requires a larger normalization constant. Such a behavior is expected, 
if a new component, e.g. of hadronic origin, sets in above 100 GeV, while at lower 
energies photons of electromagnetic origin dominate the spectrum. For a more de- 
tailed test of this hypothesis in the future, the differential energy spectrum at the 
high-energy end of the Fermi spectrum will be most useful. 

The most important consequence of the recent H.E.S.S. results is to significantly 
disfavor the jet scenario. In spite of the uncertainties in the normalization, the 
almost flat spectrum predicted in this model is in contradiction with the shape of 
the 7-ray spectrum observed by H.E.S.S. Moreover, the angular extension of the 
photon flux observed by H.E.S.S. is consistent with a point source at the AGN core 
and excludes thereby the radio lobes as sources. 

On the contrary, the shape of the 7-ray spectrum observed by H.E.S.S. agrees 
very well with the slope expected in the core model. The Fermi measurements 
restrict additionally the source model, excluding the broken-power law case that 
leads to an excessive photon fluxes in the GeV range. In summary, a primary proton 
energy dN/dE oc with a ~ 2 or smaller in the case of acceleration close to the 
core is consistent with the H.E.S.S. and Fermi observations, assuming that in the 
GeV range the photon flux has dominantly an electromagnetic origin. 

Finally, we comment on the neutrino fluxes to be expected from Cen A. Calcu- 
lating the expected event number in a neutrino telescope requires a deflnite choice 
of the experiment. Centaurus A is from the location of Icecube only visible from 
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Fig. 1. Particle fluxes from Cen A normalized to the PAO results, see the text for a description. 



above, and thus the background of atmospheric niuons allows only the use of con- 
tained events that carry essentially no directional information. By contrast, a neu- 
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trino telescope in the Mediterranean could make use of the muon signal and the 
directional information. The resulting event numbers both for cascade and shower 
event number per year observation time are summarized in Table [1] For a discussion 
on the connection between Cen A as source of UHECRs and the associated diffuse 
neutrino flux see Ref. [TOl 



Table 1. Number of neutrino events expected per year. 
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4. Conclusions 

In summary, the combination of data from PAO and H.E.S.S. is consistent with a 
scenario where the UHECRs observed by PAO are protons accelerated near the core. 
The VHE 7-rays detected by H.E.S.S. would then have a hadronic origin, namely, 
the interaction of the accelerated protons with the UV photons surrounding the core. 
This model is expected to be tested in the near future as the currently running PAO, 

H. E.S.S. and Fermi experiments increase their statistics. 
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